Quantum order by disorder in a semiclassical spin ice 
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We study the effects of quantum fluctuations in spin ice by considering an S > 1 quantum Heisenberg model 
with a nearest-neighbor ferromagnetic interaction J and a large non-collinear (111) easy-axis anisotropy D 
on a pyrochlore lattice. For a finite D ^ | J|, the low-energy physics is described by a (111) Ising model 
with additional second- and third-neighbor exchange couplings of order 0{J^ / D), generated by the quantum 
fluctuations arising from the transverse components of the exchange coupling. The extensive degeneracy of 
ground states in the _D — !> cx) limit is lifted, and a g = ordered state is selected via the quantum order 
by disorder mechanism, through a first-order phase transition at low temperatures. We propose that quantum 
dynamics in spin ice can be tuned by engineering the local single-ion anisotropy. 

PACS numbers: 75.10.Jm, 75.5().Dd,75.3().Gw,75.30.Kz 
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Introduction. — Spin ice materials are magnets W\\h ferro- 
magnetic interactions on the geometrical frustrated pyrochlore 
lattice of corner-sharing tetrahedra with strong single-ion 
anisotropy'- . The presence of local, non-collinear easy axes 
in these materials leads to pseudo Ising spins with effective 
antiferromagnetic interactions"^, and the system is highly frus- 
trated. At low temperatures, the system enters a cooperative 
paramagnetic state, and the magnetic moments(spins) obey 
the "ice rules", with two spins pointing into the center of each 
tetrahedron and two spins out. This local organizing rule gives 
rise to an exponentially large number of degenerate ground 
states, resulting in a nonzero residual entropy, and the system 
remains magnetically disordered. When a magnetic field is 
applied along the [100] direction, the spin ice enters a g = 
ordered state with saturated magnetization through a topolog- 
ical 3D Kasteleyn transition at low temperatures^'^, and the 
low-energy magnetic excitation involves a collection of spins 
on a string spanning the entire system^*^. 

An exciting direction of research in frustrated magnetism 
is to understand how the extensive classical ground state de- 
generacy is lifted by quantum fluctuations and what types of 
exotic phases may emerge. The local (111) anisotropy and 
the structure of the pyrochlore lattice in spin ice, however, 
precludes the introduction of quantum dynamics by a global 
magnetic field transverse to the Ising spins, although the off- 
diagonal terms in the dipolar Hamiltonian may suffice to in- 
troduce quantum dynamics in the rare-earth based spin ice 
materials. This issue has been addressed by either adding 
multiple-spin interactions^, or taking into account the crystal 
field excitations'". Here we propose a new route to explore the 
semiclassical spin ice by considering the quantum dynamics 
generated at finite single-ion anisotropy. 

Classical spin ice models based on Ising-Iike spins along 
the local (111) easy axisiii^^ have been very successful in 
explaining the properties of the rare-earth based spin ice ma- 
terials, such as Ho2Ti207 and Dy2Ti207. In these materials, 
the anisotropy gap A « 200 ^ SOOK^*, is much larger than 
the nearest-neighbor exchange \J\ « IK'-^, and the Ising-like 
models are well justified. However, one might expect in yet- 
to-be-discovered spin ice materials based on transition metal 



ions, where the exchange is larger and spin values are smaller, 
transverse quantum fluctuations due to finite anisotropy be- 
come important. Intriguing new ordered states may emerge 
due to the lifting of the macroscopic degeneracy via the "or- 
der by disorder" mechanismil. In the spin ice, the non- 
collinearity of the local (111) easy axes generates non-trivial 
quantum corrections at the order of 0{J^/D), in the form of 
second- and third-neighbor exchange couplings. This should 
be contrasted with the S > 3/2 triangular and kagome lattice 
antiferromagnets with a global single-ion anisotropy, where 
the leading non-trivial quantum corrections are multiple-spin 
interactions of a higher order ©( J'^/Z)^)'^"'^. 

In this paper, we consider an 5 > 1 ferromagnetic Heisen- 
berg model with a finite anisotropy D ^ \J\ along the local 
(111) easy axis. The low-energy physics is described by a 
(111) Ising model with additional second- and third-neighbor 
exchange couplings. We find at low temperature, the trans- 
verse quantum fluctuations due to finite anisotropy act as an 
order by disorder mechanism, and a q = ordered state is 
selected (Fig. [T]). This state is reminiscent of the g = or- 
dered state in Ref.^, where a 3D Kasteleyn transition is pro- 
posed. However, the specific heat in our model diverges as 
a power law, instead of logarithmically when the tempera- 
ture approaches Tc from above, and the topological spanning 
string excitations are no longer degenerate. We argue that this 
opens up new opportunities to tune the quantum effects in spin 
ice by engineering the local single-ion anisotropy, and the re- 
sulting semiclassical spin ice will provide a new playground 
to study the quantum order by disorder phenomena. 

Effective Hamiltonian. — We start with the quantum spin 
Hamiltonian (S > 1) for the nearest-neighbor ferromag- 
netic Heisenberg model (J < 0)with a local (111) single-ion 
anisotropy. 
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where i, j are the FCC lattice points, a, b label the sublattices 
inside a single tetrahedron, and the summation is over all the 
nearest-neighbor pairs, e^ corresponds to the local (111) easy 
axis for sublattice a, and D > Q confines the spins to the 



local easy axis. For D — 0, the model becomes a ferromag- 
netic Heisenberg model with a trivial ferromagnetic ground 
state. In the limit of D — > oo, the Heisenberg model maps 
into a (111) nearest-neighbor Ising model. The ground states 
obey the 2-in-2-out ice rules for each tetrahedron with an ex- 
tensive degeneracy. In the case of finite anisotropy, we ex- 
pect corrections to the classical (111) nearest-neighbor Ising 
model. In a classical ferromagnetic Heisenberg model with 
finite anisotropy, a magnetically ordered ground state with a 
four-sublattice structure is found'**. We adopt the local basis 
for each sublattice, and the local easy axis is chosen as the 
local quantization axis for the spin operators^^. The Hamilto- 
nian can be rewritten asH — Ha + Hi, where 
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Here Uab, Cab, and dab are the geometrical prefactors in the lo- 
cal basis expansion^". We construct the effective Hamiltonian 
in the finite D ^ \J\ limit following the standard degenerate 
perturbation theory^', and Hcff = PHP + PHRHP, with 
P^Ec l*o)($o| and R - E/i^o l^pK^o ^ E^)-^^p\, 
where Ea = {^a\Ho\^a)- Here (I'I'q)} are eigenstates of 
Ho- In the large D limit, the ground states of Ho contain two 
maximum spin configurations | ± 5*) at each site. 

In contrast to the 2D models with a strong global easy- 
axis anisotropy, the non-collinear local (111) easy axes intro- 
duce extra terms previously absent in the Hamiltonian, such 
as 5+5+ Si Si, S+Si Si S^ etc. For 5=1, the effective 
Hamiltonian will include transverse terms since it is possible 
to bring a spin state from | ± 1) to | =F 1) via two spin raising or 
lowering operations. For 5 > 1, it merely consists of logitu- 
dinal terms since the net effect of the perturbation cannot alter 
the spin states | ± 5) at each site. In the following we will re- 
strict our discussion to 5 > 1. In the second order term of the 




FIG. 1 . (Color online) Left: Pyrochlore lattice showing g = mag- 
netically ordered spin structure (black narrows) and a string defect 
(red arrows). Second- and third-neighbor pairs are indicated. Right: 
The energy distribution of the spanning string excitations from the 
simulation of L = 64, starting from a g = ordered state and 10^ 
strings are generated by random walks. 



effective Hamiltonian, 5+ and 5~ should appear in pairs on 
each site for non-vanishing matrix elements. Thus, nontrivial 
second- and third-neighbor interactions result from combina- 
tions of S^^S^^ and S^^Sj^ in Eq. ([T]l- Such perturbation gen- 
erates an effective Sf^S^b exchange interaction with a virtual 
spin raising and lowering process at an intermediate site k. 
All other second order terms merely renormalize the nearest- 
neighbor exchange. Therefore, extra exchange interaction in 
the effective Hamiltonian are couplings between sites linked 
by two intermediate nearest-neighbor bonds (Fig.lTJ. 

Ising equivalents. — We can map the Heisenberg operators 
into an effective Ising Hamiltonian. There are four operator 
equivalents?!!: (1) 5^ = Saf,, (2) S^Sf, = 25x (l±afJ/2, 
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where a^ = ±1. Using these operator equivalents, the new 
effective (111) Ising model is written as 

Hes = const. + Ji 



J; 
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where {{i,j)) and {{i,j)) are the second and third neigh- 
bor pairs. The Ji and J2 both are antiferromagnetic, with 

- ^"l > for 5 > 1 and 
is ferromagnetic, with J3 = 



Ji = - 
J2 
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For a given site, there are six Ji bonds, twelve J2 bonds, 
and six J3 bonds with an intermediate site (Fig. [U. When 
_D — > 00, the model reduces to the (111) nearest-neighbor 
spin ice model. The appearance of J2 and J3 is due to non- 
trivial quantum fluctuations, and we analyze their effects on 
the highly degenerate spin ice states via both mean field the- 
ory and classical Monte Carlo simulations. In the following, 
we choose a representative value of J2/J1 — 0.005, which 
corresponds to \J/D\ w 0.015 in the large 5 limit, and all the 
energies are in units of Ji . 

Mean-field theory. — To study possible ordering, we 
first analyze the effective Hamiltonian Eq. (O by the 
mean-field theory22i2i. The free energy up to the 
quadratic order at temperature T is given by J^ = 
lE,Eab{^J"'i<l)+TSa,b)m-m'Lg + 0(m4), where 
to" = jj J2i e*'*'"' (Jia, and a, b are sublattice indices. The 
matrix J'"'' (q) is the Fourier transform of the interaction ma- 
trix J'.°j^ — J{y1 — r^). The lowest eigenvalue of J^ is asso- 
ciated with the first ordered mode of the model. In the infinite 
anisotropy limit, the lowest eigenvalues form a q-independent 
flat band, which corresponds to the macroscopically degener- 
ate ground states and no preferred ordering is selected. This 
degeneracy is lifted in the presence of J2 and J3, and the low- 
est eigenvalue minimum is located at q = (0, 0, 0)(Fig. |2|, 
and the system develops a g = long-range order. 

Monte Carlo Simulation. — With the information of a pos- 
sible q — Q order from the mean-field analysis, we study 
the effective Hamiltonian Eq. ([3]) using Monte Carlo simu- 
lations. The simulation is done on the pyrochlore lattice with 
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FIG. 2. (Color online) Lowest eigenvalue of the interaction matrix 
e7"''(q) at wave vectors q = 2-K{h, h, I). The minimum occurs at 
q= (0,0,0). 
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FIG. 3. The energy histogram for L = 8 near the transition temper- 
ature. J2 — O.OO5J1 and J3 — — 2J2. The double peak structure in 
the histogram indicates a first order phase transition. 



periodic boundary conditions measuring L cubic unit cells in 
each direction, which amounts to a total of iV = 16 x L^ 
spins. We perform our simulations using the parallel tem- 
pering algorithm^^, so that the simulation can reach ergodic- 
ity more efficiently at low temperatures. The simulations are 
conducted in parallel at a series of temperatures, and swaps 
of configurations between these parallel simulations are pro- 
posed, allowing the low temperature system of interest to es- 
cape from local free energy minima where it might otherwise 
be trapped. In addition, the configurations still favors ice rule 
for Ji ^ J2 in each serial run, so we also employ loop up- 
dates to avoid the ice-rule breaking energy barrier^*. In our 
simulations, for L = 4 to 7, we carry out 2000 configuration 
swaps and 200 loop updates between swaps in both equilib- 
rium and sampling processes; for L = 8 to 10, 4000 con- 
figuration swaps and 200 loop updates between swaps during 
equilibration and 2000 swaps and 200 loop updates between 
swaps in sampling are carried out. Figure [3] shows the en- 
ergy histogram for L = 8 at temperatures near the transi- 
tion. The double peak feature at the transition indicates the 
coexistence of two distinct phases at Tc. Figure |4] shows the 
temperature variation of the specific heat for various lattice 
sizes. As the size increases, the peaks grow and peak widths 
decrease. There is also a shift in the temperature correspond- 
ing to the specific heat peak Tp. The specific heat in a first 
order transition should be extensive, and diverges in the ther- 



FIG. 4. (Color online) Variation of specific heat as a function of tem- 
perature for several different sizes. Insets: (a) Specific peak versus 
L^. A clear extensive behavior is shown in the linear fit. (b) The 
peak temperature versus 1/L^. The transition temperature in ther- 
modynamic limit is extrapolated to be T^ = 0.181. The behavior 
is consistent with the properties of a first-order phase transition. 
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FIG. 5. (Color online) Scaling of the specific-heat data. Tp corre- 
sponds to the temperature of the specific heat maximum for a given 
size L. 



modynamic limit. We fit the size dependence of the specific 
heat by C == 2.142.3 ^ o.l26 (inset (a) of Fig. gli. On the 
other hand, the peak temperature T^ for a given size L can be 
fit to T^ = T^ + 1.78/L3 with T^ = 0.181 (inset (b) of 
Fig. nil. The above two relations suggest that it is possible to 
perform scaling analysis on the specific heat data using C and 
AT = T — Tf^. Figure|5]shows the scaling of the specific heat 
data, and the scaling works reasonably well. We conclude that 
this transition to the q — Q ordered state is a first order phase 
transition^^. 

At low temperatures, the ground state corresponds to a g = 
ordered state and carries saturated magnetization toward one 
of the [100] directions(Fig.[Tll. The ground state energy can be 
computed exactly, Eg/N = — Ji — 2, J 2 + 3 J3 = — Ji — 8 J2. 
In the ice-rule states, the low energy excitation is a collection 
of flipped spins on a string, as a single spin-flip costs higher 
exchange energy of order 0(Ji). The spin ice state and the 
(7 = state both obey the local 2-in-2-out constraint but their 
excitations show distinct topology. In the spin ice state, the 



strings form loops of finite lengths, while in the q — ordered 
state, the string will span the entire system. These spanning 
strings excitations carry energy and entropy proportional to 
the segment length of the string. At high enough temperature, 
such excitations will lower the free energy and the ordered 
state is destroyed. This is similar to the case of the spin ice in 
a [100] field^, where a 3D Kasteleyn transition into a g = 
ordered state at low temperature is proposed. In our model, 
however, the specific heat diverges as a power law when the 
temperature approaches T^, from above; while in the case of 
a spin ice in [100] field, it diverges logarithmically. We note 
that in our model the exact degeneracy of the spanning string 
excitations in Ref.^ is lifted and the excitation energy depends 
on the path the string traverses. Figure [T] shows the segment 
energy distribution of the spanning string excitations from the 
simulation of L — 64, starting from a q — ordered state 
and a total of 10^ strings are generated by random walks^. A 
broad distribution is clearly observed due to the path depen- 
dence of the energy in the string creation. 

Conclusion. — We propose that a semiclassical spin ice 
with finite anisotropy provides a new playground to study 
the quantum order by disorder phenomena. Transverse quan- 
tum fluctuations, which can be tuned by engineering the lo- 



cal single-ion anisotropy, generate additional second- and 
third-neighbor exchange interactions of order 0{J^ /D) in the 
(111) Ising model. These new interactions lift the extensive 
ground state degeneracy in the infinite anisotropy limit and 
select a six-fold q ~ ordered ground state carrying satu- 
rated magnetization toward one of the [100] directions. Al- 
though the topological characteristics of the string excitations 
are similar to those in Ref.^, the critical behavior is quite dif- 
ferent due to the path dependence of the spanning string en- 
ergy. Interesting questions remain on the effects of an external 
magnetic field and dilution^*^, which requires further study. 
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